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SUMMARY 


A nonlinear least squares algorithm for aircraft parameter 
estimation from flight data was developed. The postulated model for the 
analysis represented longitudinal, short period motion of an aircraft. The 
corresponding aerodynamic model equations included indicial functions 
(unsteady terms) and conventional stability and control derivatives. The 
indicial functions were modelled as simple exponential functions. The 
estimation procedure was applied in five examples. Four of the examples 
used simulated and flight data from small amplitude maneuvers of the 
F-18 HARV and X-31A aircraft. In the fifth example a rapid, large 
amplitude maneuver of the X-31 drop model was analyzed. From data 
analysis of small amplitude maneuvers it was found that the model with 
conventional stability and control derivatives was adequate. Also, 
parameter estimation from a rapid, large amplitude maneuver did not 
reveal any noticeable presence of unsteady aerodynamics. 
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SYMBOLS 


Q* fljji Clq 

a x’ a z 
B(X’ Bq 

1>1 

^L’ Qn> 


C«, Cq 

C La (t),C Za (t),CmJt) 
Cov( 0) 

C y Cft, Cq 

* * 

C<x> Cq 

C 

F(t),F L Jt) 

g 

h-ho 

H 

J x < Iy< h 

hz 

J(6) 

£ 

k 

m 


parameters in indicial functions 

longitudinal and vertical accelerations, g units 
parameters defined in table I 

parameter in indicial function, Usee 
lift, pitching-moment and vertical-force 

coefficients 

parameters defined in table I 
indicial functions 

parameter covariance matrix 

parameters in indicial functions 

parameters in indicial functions defined in eq. (23) 

mean aerodynamic chord, m 
deficiency functions 

acceleration due to gravity, m/ sec 2 
location of aircraft center of gravity from the 
aerodynamic center of the wing expressed in c 
sensitivity matrix 

moments of inertia about longitudinal, lateral, 

and vertical body axes, kg-m 2 
product of inertia, kg - m 
cost function 
characteristic length, m 
tail arm, m 

mass, kg 
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number of data points 
number of parameters 

roll rate, pitch rate, and yaw rate, rad /sec or 
deg /sec 

dynamic pressure, pV 2 /2, Pa 

coefficient of determination 

wing area, m 2 

standard error 

time lag, f/V, sec 

time, sec 

airspeed, m / sec 

state variable in eq. (10) 

dependent variable 

measured dependent variable 

angle of attack, rad or deg 

control surface deflection, rad or deg 

canard and thrust vectoring vane deflection, rad 

or deg 

vector of unknown parameters 
unknown parameter 

a) measurement error 

b) error due to approximation (Appendix A) 
air density, kg/ m 3 

time delay, sec 

angular frequency, rad /sec 



Superscripts: 


Subscripts: 

t 

w 

a.c. 


derivative with respect to time 
estimate 

tail 

wing 

aerodynamic center 


Matrix exponents: 
T 
-1 


transpose matrix 
inverse matrix 


Derivatives of aerodynamic coefficients C a where the index a = L 
_ dC a _ dC a _ dC a _ dC a 

a “~ d V_ a «- da’ as 95’ d & 

2V 2V 


derivatives and C* defined in eq. (22) 

a a U S 

derivatives L a ^ q § and M aaq § defined in table I. 


, Z, or m 
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INTRODUCTION 


Possible effects of unsteady aerodynamics on aircraft motion were 
investigated as early as in the forties and fifties (see reference 1 and 2). So 
far, however, there have been only a few attempts to analyze flight data with 
the intent to estimate parameters characterizing unsteady aerodynamics. 
In reference 3, a procedure for estimation of unsteady aerodynamic terms 
was proposed and applied to longitudinal data of a scaled model. In this 
procedure, the conventional stability and control derivatives were estimated 
first. Then, the resulting residuals were used for estimation of unsteady 
terms. Reference 4 addressed the formulation of indicial functions in the 
equations of motion and the identifiability of these functions. Examples 
were given using simulated data and then flight data of a drop model from 
the post-stall and spin region. It was found that the unsteady part of the 
model was not identifiable for the responses below the angle of attack of 40°. 
In reference 5 a simple vortex system was used to incorporate the unsteady 
aerodynamics into the longitudinal equations of motion. Neither simulated 
nor flight data of a general aviation aircraft undergoing small amplitude 
maneuvers at low angles of attack showed any significant difference in 
parameter estimates using either the model with or without unsteady 
terms. Finally, in reference 6 and 7 an unsteady aerodynamic model was 
determined from longitudinal, large amplitude maneuvers. In this case 
the modeling was based on internal state variables rather then indicial 
functions. From the investigation, however, it was not clear how the model 
with time-invariant parameters would explain the measured response. 

The present report is the third part of the series devoted to the 
modeling of unsteady aerodynamics in the equations of motion and to 
estimation of aerodynamic parameters from experimental data. In the 
first part of this series (reference 8), linear models for aerodynamic forces 
and moments were formulated. In the second part (reference 9), developed 
models were used in the analysis of wind tunnel data from forced 
oscillation test. The purpose of this report is to postulate an aircraft model 
with a simple form for the unsteady aerodynamics, investigate the effect of 
unsteady terms on the motion of an aircraft and estimate the aerodynamic 
parameters in the postulated model from flight data. Because the report is 
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considered as a preliminary study into the problem, all examples 
presented, except one, are limited to the longitudinal, short period motion 
with linear aerodynamic equations. Only the last example investigates the 
adequacy of a model without unsteady terms for a longitudinal, large 
amplitude maneuver. The report starts with the formulation of model 
equations followed by the development of a parameter estimation algorithm 
based on the least square procedure. Then, the algorithm is applied to 
simulated and flight data. The results obtained are discussed and 
concluding remarks drawn. 
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EQUATIONS OF MOTION 


As an example of aircraft motion with unsteady aerodynamic terms 
the short period longitudinal motion will be studied. The development of 
these equations is similar to that in reference 8 . The short period 
longitudinal motion can be described as 


a = q-£™C L {oc(t),q(t),S(t)) 
2m 

2 — 

q = ^ 7 —C m (a(t),q(t),8(t)) 

Zly 


It is further assumed that the aerodynamic model equations can be 
formulated as 

C L (t) = \C Lq (t-r)-^a(r)dT + ^C Lq (°° )q(t) + C Lg (°° )S(t) 

° <5 

C m ( t ) = J C ma (t - r)-^ a( r)dr + ^C mq (« )q{t) + C mg (<*> )8(t) 

o 

where C L Jt) and C m Jt ) are indicial functions representing the response 
in C L and C m to a unit step in a, and the remaining coefficients are the 
conventional stability and control derivatives. In a more general case they 
might be replaced by indicial functions. 

The indicial functions in equation (2) achieve steady state values 
C L (°o) and C ma (°°) respectively. Because of this property, the indicial 

function C L (t) can be expressed as 

C ta O) = C io (~)-F ln (i) G 

where F La (t ) is called the deficiency function (ref. 2). A similar 
relationship holds for C ma ( t) . 

The most difficult part of modeling aerodynamics in eq. (1) is the 
postulation of expressions for indicial functions. If the model is to be used 
in parameter estimation from experimental data, then it should be 
parsimonious. This means that the model must fit the data well with the 
smallest number of parameters. The problem of formulating aircraft 
indicial functions was addressed in references 1, 2, 5, 10 and 11 . For this 


7 



study, the approach of references 10 and 11 was used. Reference 10 extends 
the two-dimensional potential theory of an airfoil in nonuniform motion to a 
wing with finite aspect ratio and elliptic spanwise loading. Further 

extension to tapered, swept wings in incompressible flow is covered in 
reference 11. It is shown that in all cases the indicial function C^ a (t) can 

be formulated as 


C La (t) = a( l-e _6l *) + c 
= CL a (°°) — ae blt 


(4) 


which was also the form used in ref. 9 for the analysis of wind tunnel data. 

The indicial function for a wing-tail combination, however, becomes 
more complicated. The indicial function C^ a (t) must be considered to 

represent the combined effect between the wing and tail. In linear theory 
the resulting indicial function is equal to the sum of the components (see 
ref. 2): 

(1) the indicial function of the wing alone, the tail being at zero angle of 
attack; 

(2) the indicial function of the tail alone, the wing being at zero angle of 
attack; 

(3) the response in lift of the tail to a step change in angle of attack of the 
wing, the tail being at zero angle of attack; 

(4) the response in lift of the wing to a step change in angle of attack of 
the tail, the wing being at zero angle of attack. 

The last two components express the interference effects. 

Reference 2 presents a detailed qualitative analysis for supersonic 
speed. References 1 and 5 cover the subsonic case by concentrating on the 
effect of the wing wake on the lift of the tail and neglecting the influence of 
the tail on the wing. Examples from reference 10 and 1 are reproduced in 
figures la and lb. In these figures the indicial function of an elliptic wing 
with aspect ratio 6 and the indicial function of a tail with aspect ratio 3 are 
plotted against the nondimensional time, Vt/(c/2) . As seen on figure lb 
the lift on the tail becomes infinite when the leading edge of the tail reaches 
the position of the starting vortex. After that, the lift decreases and, with 
increasing time, approaches a steady value. Although the indicial lift 
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shows an infinite value, the convolution with the angle of attack results in 

finite lift at all points. The dotted line in figure lb is an approximation to 
C L( (t) representing the development of the lift on the tail by a single lag 

function. However, the lag is greater than that indicated by the tail length. 

The resulting indicial functions for the lift and pitching moment of 
the whole aircraft can be obtained from the expressions 

Cl<t) = C K (t)+^C Lt (t) (5) 

C„(O = C m "-Hh-h 0 )C K (t)-t&-C L ,(O (6) 

In more detailed analysis, should be replaced by the lift of a wing-body 

combination. Using the results in figure la and lb, the indicial function 
C m (t) and its approximation were computed from equation (6) for 

r = o, h - h n = 0.1 and the tail volume equal to 0.6. Both functions are 

m a.c U 

plotted in figure lc where the dotted line indicates the approximation. 

For the following analysis the indicial functions are approximated by 
rather simple forms as 


C L Jt) = C L J°°)-a a e b * 

C„ a (0 = C„ o (°°) - a q e~ 1 ’' 1 

which means that the pitching-moment indicial function may be suitable 
for a flying wing and perhaps for an aircraft with a short tail length. 
Substituting eq. (7) into eq. (2) and integrating by parts results in 
aerodynamic equations 


C L (t) = c a a{t) + a a b } J e b ‘ r a(t- f )dr + C L q(t) + C Lg 8(t ) 

0 v 

C m (t) = Cq(x(t) + ciqbjje 1 oc(t — T)dr+ — C m ^q(t) + C m ^8(t) (8) 
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where the simplified notation for aerodynamic derivatives was used. 
Combining eq. (1) and (8), and introducing dimensional parameters, the 
equations of motion can be written as 

t 

a = -C a a(t) - B a J e~^ r a( t - r)dr+(l -L q )q- L s 5 

0 

t 

q = C q a(t) + B q J c _ ^ T af t - r)dr + M q q + M s 8 

o 


( 9 ) 


where the parameters in these equation are defined in table 1. Introducing 
a new state variable 

x a = je~ b,r a(t-T)dT 
0 


and the corresponding state equation for this variable 

x a = a-bi x a 

equations (9) can be rewritten in the state-space form as 


a 


1 

1 

sP 

1 

i 

** 
i 

a 


\ 

i 

i 

q 


Cq Mg Bg 

q 

+ 

M S 

Xa. 


1 

1 

o 
1 

_ x a_ 


0 
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If no unsteady aerodynamics are considered, eq. (2) are simplified as 

C L ( f ) = C L a «( 0 + ^ ( C Lg ®(0 + c L q ^ )) + C L S 5 (0 

C m (0 = C ma cc(t) + —^C m ^dc(t) + C m ^q(t)^ + C ms 8(t) (11) 


and the longitudinal short period equations will take the form 


1+L a 0 

-M a 1 


La 1 L q 


-L S 

M S 


( 12 ) 


'a "+q 

where the parameters L a , M a , L a and M a are also defined in table 1. 

Equations (11) represent an approximation to the more general form 
given by eq. (2). The difference between those two formulations can be 
interpreted as an error due to omitting the unsteady aerodynamic terms. 
The dependence of this error on the aircraft geometry and motion variables 
is investigated in reference 4 and Appendix A. It was found that for a 
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is investigated in reference 4 and Appendix A. It was found that for a 

harmonic motion with given frequency (0, the relative error is equal or less 

than the reduced frequency of the motion, i.e. 

co£ 

relative error < — 

V 

where l is a characteristic length. This result indicates that the error will 
increase with increased frequency of motion and characteristic length and 
decrease with the airspeed. During transient longitudinal motion, the 
prevailing frequency is close to the natural frequency of the aircraft. The 
characteristic length will depend on the physics of the flow which creates 
unsteady effects. For unsteady lift on the wing l-c , for the interaction of 
the wing and the tail i is equal to the tail length, l = £ x . When fuselage 

vortices are interacting with the tail, £ will be equal to the fuselage length. 
However, the airspeed during the maneuver will have the main effect on 
reduced frequency. Therefore the unsteady effects will be more pronounced 
at low airspeed which means at high angles of attack. 

PARAMETER ESTIMATION ALGORITHM 


A regression method for estimation of aerodynamic parameters from 
measured input/output time histories was selected because of its simplicity. 
Considering aerodynamic model equations (8), the regression equation for 
each of the model equations can be written as 

t _ (j £ 

z(i) = 6 ] a(i) + d 2 d 3 \e 3t a(ti - T)dr+ 0 4 —q(i) + 0 5 S(i) + e(i) 

o v v ' 

i = 1, 2, N 

In this equation 6j are the unknown parameters, z(i) is the measured 
independent variable, e(i ) is the measurement error and N is the number 
of data points. For the pitching-moment coefficient 

— c q i @3 = b] , @5 = Cm s 

@2 = a q > ~ Crriq 

and z(i) = C m (i). The parameter estimates are obtained by minimizing the 


cost function 


J= I £i 2 

i=l 


N 2 

= I [z(i)-y(i)] 

i=7 


where y(i) is the computed independent variable 


( 14 ) 
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Because regression equation (13) is nonlinear in the parameters the 
minimum of the cost functions is computed from the linearized form 


N 

J= I 

i=; 


n p dv( i ) 

z(i)-y 0 (i)- 2 


i2 


j=l Mj 


(15) 


where dy/ddj are the sensitivities 


l 

= b l je~ b,r a(t - rjdr, 


*L = o, J* 

902 u 

= a q j(\ - b l r)e~ b]T a( I - T )dx . 

3y 


0 

= 1 
dd 4 ~ v q 


dy_ 

dQs. 


= d, 


The sensitivities and the independent variable yg are computed for nominal 
values of the parameters. Minimization of (15) leads to an iterative 
procedure where the estimates after the rth iteration are equal to 

0 T+] =e T + A0 T (16) 

The parameter updates A0j are obtained as 


A0 = 


N „ 

I#/#, 

i=l 


-l -1 


I Hj{z(i)-yo(i)) 
J « =1 


(17) 


where H is the sensitivity matrix 

dy 


H = 


dOj \ ’ 


j = 1,2,. 


„n t 


Assuming that a. , q , and 8 are measured without errors and that e 
forms a random, white sequence, the parameter covariance matrix can be 
estimated as 


Cov(<9) = , 


N „ 

I HjHi 
i=l 


-1-7 


(18) 


where 


, 2 = 


_ '(®) 


N -n r 


For the aerodynamic equation expressed as 


d f 

C m (t) = C ma a(t)-a q \e- b ' T —cc(t-T)dT + -C mi/ q(t) + C m5 5(t) 


(19) 
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( 20 ) 


the corresponding regression equation is 

z(i) = Oja(i) - 0 2 j e~ 03T -4- a( t t - r)dr + 04 ^-q( i) + #5<5 ( 0 
o dr V 

and the sensitivities and are 

-^- = -\e~ b[T — a(t- r)dr 
de 2 J dr 

-^- = a a f r e~^ x -a(t- r)dr 
de 3 q l dr 

If the parameter, 0 3 , in equation (13) or (20) is assumed to be known, the 
estimation of the remaining parameters is reduced to an ordinary least 
squares method. 

A different estimation scheme based on the least-squares principle 
was introduced in reference 12. It transforms state equation (1) and 
equation (2) into the complex domain. Then the modulating function 
technique of reference 13 is applied to those equations in order to obtain a 
computationally suitable form of the regression equations. These equations 
are linear in the parameters. The estimates of aerodynamic parameters in 
(2), however, require a solution of nonlinear algebraic equations. These 
equations could be rather complicated even for a small number of unknown 
parameters. For that reason the algorithm of reference 12 was not used in 
this report. 


ANALYSIS OF SIMULATED DATA 

Two examples of simulated data were generated to explore the effects 
of unsteady aerodynamics on system output variables and aerodynamic 
parameters estimated from the input/output data. Both examples are of 
aircraft performing small-amplitude longitudinal maneuvers from trim 
conditions at 30° angle of attack. The model used to generate the data was 
that given by equations (10) and (12). In addition to the computed output 
variables, a and q, the time histories of the lift and pitching moment 
coefficients with unsteady aerodynamic terms were also computed from eq. 
(8). Three different approaches were used for estimation of unknown 
parameters. In the first approach the aerodynamic model with only 
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conventional stability and control derivatives (see eg. (11)) was considered. 
In the remaining two cases, the correct model given by equation (8) was 
applied with the parameter bj either fixed or estimated. Thus the first and 
second approach resulted in a simple linear regression. For the third 
approach, the estimation algorithm described by equations (15) to (18) was 
applied. 

Advanced Fighter Aircraft 

Aircraft characteristics and flight conditions for an advanced fighter 
aircraft are summarized in table II. These characteristics are close to 
those of the F-18 High-Angle-of-Attack Research Vehicle (HARV) as 
reported in reference 14. The value of parameters c, a and bj in the 

indicial functions are included in table III and IV. The computed time 
histories of the angle of attack, pitch rate and two coefficients C L and C m 

are plotted in figures 2 and 3. Figure 2 also includes the input variable S(t). 
Both the response variables and aerodynamic coefficients show only small 
differences between their values from equations with and without unsteady 
aerodynamics. 

For parameter estimation, measurement noise was added to the 
coefficient Q, and C m computed with unsteady terms. The noise was white 

and Gaussian with zero mean and variance g 2 =6.9x10 3 for and 
1.55X10 -4 for C m . The results from three difference cases are given in tables 
III and IV. Included are parameter values and their standard errors, 
coefficients of determinations, R 2 , which define the amount of information 
in the data explained by the model, and the standard errors of the 
coefficients, s(C L ) or s(C m ), which represent fit errors. The estimation 
algorithm which includes bj as an unknown parameter did not converge 
where the data with the measurement noise were used. Convergence was 
obtained only with the reduced amount of noise. In the given example the 
noise was reduced to 1/20 of the original measurement noise in C L and to 
1/3 of the original value for C m . 

The model with conventional stability and control derivatives fitted 

the data well and the estimated parameters were close to their true values. 
One should remember that, in this case, the estimate of C m represents the 
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sum of C m and C m with the true value of -12.5. The remaining estimates 

rn q ni oi 

of C m and CL are also affected by terms containing C m . . Similarly, the 
parameter estimates of C^ a , and C^ 5 contain a small contribution due 
to C L . When the unsteady model was used with the parameter b\ fixed, no 

noticeable improvements in R 2 , s(Ci) or s(C m ) were observed. There was, 

however, an increase in the standard errors of parameters accompanied by 
high correlation between parameters c a and C Lq , a a and c a , and a a and 

Ci as can be seen in table V. The same applies for the correlation of 

parameters in the pitching-moment equation. By adding the parameter b\ 

to the unknown parameters, the correlation of parameters become even 

worse than in the previous case. As indicated by table VI, all three 
parameters c a , a a and b\ , or Cg, a q and b\ are highly correlated. 

Improvement in the accuracy of the estimated parameters is a result of 
decreased measurement noise variance. 

Experimental Aircraft 

The aircraft characteristics and flight conditions for this example 
are presented in Table VII. The characteristics were selected close to those 
of the X-31A aircraft. This aircraft has three longitudinal controllers, the 
tailing edge flaps, canard and thrust vectoring. In the model developed, it 
was assumed that the control system generates the canard deflection, 8 C , 
proportional to a, and thrust deflection, 8^, proportional to tailing edge 
flaps deflection, 8. The lift coefficient expressed with conventional stability 
and control derivatives has the form 

C L (t) = C La a(t) + ^(c L . a(t) + C Lq q(t)) + C Ls S(t) + 8 c (t) 

+ C L S &t vW 

d tV 

= c' La a(t) + ^(c L .d(t)+C Li q<tj)+ c' Ls S(t) 

where 

C L a = C L a +k <x C L Sc 

C*L S = c Ls + k s C L ^ 


( 21 ) 


( 22 ) 
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Similar expressions hold for C m (t). The aerodynamic equations with 
unsteady aerodynamics were obtained from equation (8) with c a and c q 


replaced by 


-a 


= Ci -a, 


a 


( 23 ) 


a. 


"q ~ w m a 

The values of parameters a a , a q and b \ , were estimate from wind tunnel 


oscillatory data as described in reference 9. The values of these parameters 
are in the second column of tables VIII and IX. 

As follows from figures 4 and 5, the effect of unsteady aerodynamics 
is now more pronounced than in the previous example. For parameter 
estimation, the time histories C L (t) and C m (t) computed from modified 


equations (8) were again corrupted by measurement noise with the 
variance of 4.0xlCT 3 for C L and 1.19XKT 4 for C m . The increased 
significance of unsteady aerodynamics was reflected in the accuracy of the 
estimated parameters and in the value of R . When compared with 
estimates related to equation (21). The estimation with the correct model 
and fixed parameter b\ moved the parameter mean values closer to the true 
values, but there was an increase in parameter standard errors. There 
was some improvement in R 2 in the lift equation, whereas the coefficient of 
determination in the pitching-moment equation changed only slightly. 
Finally, where all parameters were estimated their accuracy decreased 
and high correlation between parameters of the indicial functions 
appeared. These results are all given in tables VIII and EX. 


ANALYSIS OF FLIGHT DATA 


Several longitudinal, short period maneuvers of the F-18 HARV and 
X-31A aircraft were analyzed. The purpose of this analysis was to estimate 
parameters in the model with and without unsteady aerodynamics and to 
decide whether a simple model with stability and control derivatives is 
adequate for each of these maneuvers. Reported are, however, only two 
examples, one for each aircraft. Conclusions from this set of data are 
similar to those from the remaining maneuvers. The third example, based 
on data of the X-31 drop model, was added to demonstrate the adequacy of a 
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model without unsteady terms for a longitudinal, high amplitude 
maneuver. 

Prior to parameter estimation the measured data were checked for 
their compatibility and estimated bias errors (constant offset and/or scale 
factors) were removed from the data. Then, the aerodynamic coefficients of 
the vertical force and pitching moment were computed as 
r - m S n 

c z -— a 2 


C -_2- a- 
- — n— 9 


1 9 U / r? / 2 2\ 

z — pr — p-(r - p ) 

Iy ly 


where the angular acceleration in pitch was obtained by numerical 
differentiation of measured pitching velocity. For parameter estimation the 
coefficient Cz rather than Cp was used. It was assumed that the accuracy 
of calculated Cz would be higher than that of C L because that lift coefficient 
is computed from time histories of two accelerations, angle of attack and 
thrust. 


F-18 HARV 


The measured time histories of longitudinal variables from a 
transient maneuver about a ~ 28° are given in figure 6. This maneuver 
was excited by an optimal input implemented by the pilot (see reference 15) 
and was performed without thrust vectoring. The measured coefficients 
Cz and C m are presented in figure 7. For parameter estimation, a model 

without unsteady aerodynamics and a stepwise regression of reference 15 
were applied first. An adequate model for Cz contained only two 
regressors, a and q , and resulted in coefficient of determination, 

R 2 = 97.53. For the coefficient C m an adequate model included three 

regressors, a, q and 8, and the coefficient of determination was equal to 
95.63. The predicted coefficients C z and C m from this analysis are 

compared with those measured in figure 7 where the corresponding 
residuals are also plotted. Both the values of R 2 and the residuals indicate 
that a simple model can explain a substantial part of variation in the data. 
Regardless of these findings, parameters in the unsteady model given by 
equation (8) were estimated for b i fixed on several values ranging from 0.5 
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to 1.5. No improvements in R 2 and in the residuals were observed. In 

addition, the parameter accuracy was degraded and strong multiple 
correlations among the parameters a and c, and Cz q or C m<? were 

observed. 

X-31A Aircraft 

The longitudinal maneuver presented was excited at a =68° by three 
controllers, trailing edge flaps, canard and thrust vectoring. The 
measured time histories of three output and three control variables are 
plotted in figure 8. The measured coefficients C% and C m are given in 

figure 9. Parameters were estimated in the same way as in the previous 
example. The predicted coefficients from estimated stability and control 
derivatives are compared with measured values in figure 9. The 
corresponding residuals are also included. The values of R for and 
C m were equal to 95.15 and 96.89 respectively. The residuals in C% include 

some low frequency components caused probably by external disturbances. 
Parameter estimation in the model with unsteady terms did not increase 
the R 2 value and did not significantly change the residuals. Because the 
residuals in C m resemble a random, white sequence, and because of the 

high value of R 2 , no further parameter estimations in the pitching- 
moment equation with different models was attempted. 

X-31 Drop Model 

The X-31 drop model is a 27% scale model of the actual aircraft. A 
brief description of the model and testing procedure can be found in 
reference 17. The measured data containing time histories of three output 
and two input variables are presented in figure 10. The data indicate that 
in the maneuver, the angle of attack was changed from 22° to 75° and back 
to 20° within 3.5 sec. The maximum pitch rate reached was 80°/sec which 
corresponds to 42°/sec for the full scale aircraft. The measured 
aerodynamic coefficients C%, and C m are shown in figure 11. The lift 

coefficient was obtained from measured data as 

C L = ~z^{a, x sin a-a z cos a) 
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Possible nonlinearities in the aerodynamic model equations without 
unsteady aerodynamics were modeled by polynomial splines, see reference 
18. These equations were formed as 


C a =C a (a)\ 


q=8 = 0 

S c =~ 40° 


+ c n 




a = Z, L, or m 

where 

C a (a), ... , C a$ (a) are polynomial splines of the first or second order, 
A5 C =S C + 0.698 

The stepwise regression was used to determine an adequate model for each 
coefficient from postulated spline terms and estimated parameter in these 
models. The predicted aerodynamic coefficients and residuals are plotted 
in figure 11. The values of the coefficients of determination were equal to 
99.81, 98.79 and 99.81 for Cz, C L and C m respectively. Both the plots and 

R 2 values indicate that the models with time-invariant parameters explain 
almost all variation in measured data. 


DISCUSSION OF RESULTS 


Models for the longitudinal, short period motion used in this report 
contained two indicial functions, Cl it) or Cz a (t), and C ma (t). The 

remaining terms in the aerodynamic model equations were represented by 
conventional stability and control derivatives. The analytical expression for 
the indicial function was selected in the form of a single exponential 
function with three parameters, a, c, and £q, see eq. (4). Values of the first 

two parameters are constrained by the condition that their sum is equal to 
the derivative Cz , Cl u , or C ma . A priori values for parameters a and 

can be obtained from theory or wind tunnel experiment. As indicated 
earlier, the form of indicial functions used in simulation and flight data 

analysis can approximate unsteady aerodynamics of a wing alone. The use 
of the same form for the indicial function C ma ( t ) of a wing-tail combination 

is questionable. It can be substantiated only as a first approximation before 
more rigorous analysis is attempted. 

Examples with simulated and measured data revealed that it might 
be difficult or even impossible to detect unsteady effects in data from 
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longitudinal, small amplitude maneuvers for two main reasons: these 
effects can be small enough so that they cannot be distinguished from the 
measurement noise and/or unknown modeling errors; even if there are no 
modeling errors and low measurement noise in the data, the identifiability 
of parameters in indicial functions is in doubt. Further, the examples 
showed that a model with conventional stability and control derivatives 
could be adequate. The model with unsteady aerodynamics did not improve 
the fit to the data. In addition, its use resulted in parameter estimates with 
low accuracy and high correlations. From the limited theoretical analysis 
and small amount of examples it can be concluded that, in general, 
parameter estimation might detect unsteady aerodynamic effects in 
measured data from small amplitude maneuvers if the transient motion is 
initiated at low speed and contains high frequency components, the 
characteristic length is large, and the aircraft has high gradients 
dC l / da and dC m Ida. It is also important that the measured data were 

corrected for bias errors and have a high signal-to-noise ratio. 

The last example of this report dealt with a rapid, large amplitude 
maneuver. It was expected that the effect of unsteady aerodynamics would 
be detected. The model with a small number of time-invariant parameters, 
however, explained almost all the variation in the data. This indicates that 
even that type of maneuver can not guarantee a noticeable presence of 
unsteady aerodynamics. 

CONCLUDING REMARKS 

A nonlinear least squares algorithm for aircraft parameter 
estimation from flight data was developed. The postulated model for the 
analysis represented longitudinal, short period motion of an aircraft. The 
corresponding aerodynamic model equations included indicial functions 
(unsteady terms) and conventional stability and control derivatives. In 
formulating the analytical form for the indicial functions, two conflicting 
requirements had to be taken into account: parameter estimation requires 
a simple model with a small number of parameters in order to improve or 
insure their identifiability; indicial functions should be good 
approximations to complex physical phenomena associated with unsteady 
and separated flow. The model proposed is this study was formed as a 
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simple exponential function which can approximate unsteady 
aerodynamics of a wing alone, but is questionable for the indicial function 
in the pitching-moment equation of a wing-tail combination. 

The estimation procedure was applied in four examples to simulated 
and flight data from small amplitude maneuvers of the F-18 High-Angle-of- 
Attack Research Vehicle (HARV) and X-31 A aircraft. Then, in the fifth 
example a rapid, large amplitude maneuver of the X-31 drop model was 
analyzed. From brief theoretical study and examples, the following 
conclusions can be drawn: 

1. A possibility for detecting unsteady effects in measured data will be 
enhanced if the transient motion is initiated at low speed, the 
response contains high frequency components, aircraft 
characteristic length is large, and the angle-of-attack gradient of the 
lift and pitching moment are high; 

2. from data analysis of small amplitude maneuvers it was found that 
the model with conventional stability and control derivatives was 
adequate; 

3. the model with unsteady aerodynamics did not improve the fit to the 
data over the simple model without unsteady terms. More complex 
models resulted in parameter estimates with low accuracy and high 
correlations; 

4. parameter estimation from a rapid, large amplitude maneuver also 
did not reveal any noticeable presence of unsteady aerodynamics. 
Despite the negative results from this study, more extensive research 

is needed to determine how important the modeling of unsteady 
aerodynamics is on aircraft motion. 
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APPENDIX A 


Approximation of Unsteady Aerodynamic Terms 


As an example, the pitching-moment coefficient is formulated with 
and without the unsteady term as 

T 

C m =C mo +C ma a-jF(t- x)d( r)dr +^C mq q (A 1 ) 

0 


and 

C m ~ Cjuq + Cm a a+ y {^' m a ® + ^ m q (A2) 

where the time T in the upper limit of the integral can be viewed as a delay 
after which the indicial function approaches its steady value. Thus T is 
proportional to the time the flow needs to travel a characteristic distance (, 
which means that T = £/V . 

The difference between eq. (Al) and (A2) represents an error due to 


approximation 


e'= jF(t-x)d(x)dx + ^C m6 <x 


(A3) 


It has been shown in ref. 7 that, under some simplifying assumptions, the 
counterpart to C m is proportional to the area of the deficiency function, that 


is 


c m„ 


(A4) 


substituting (A4) into (A3) 

T 


Since 
then also 


e' = ^F(t-x)a{x)dx-a^F(x)dx 

o o 

T 

F(f)[a(£- x)-a{t)]dx 


\cc(t x) a( x)\ < x\d niax (t)\ 
\d(t-x)-a(x)\< rja^f/)! 


(A5) 
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Using the last inequality the error can be further expressed as 

T 

e' <1 J|F(r)||d(£- r)-d(T)|dr 
0 

T t 

< f| F(z\i\a 

max (t)\dr<T\d max 1 |F( T)jdT 
0 0 

The maximum value of the integral in eq.(Al) is 

0 


which is used as a normalizing factor in the expression for relative 

< < l^max^l rp 

" , ... |«max(*)| 

0 

For a simple harmonic motion, the relative error takes the form 


where (o is the angular frequency 


(A6) 


error 

(A7) 


(A8) 
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Table I. - Definition of parameters in equations. 


L - pSV Cr 
a 2m Ch '* 

M - pV Sc C 

Ma ~ 2ly Cm - 

T. -P Sc Cr 
a Am Cl ° 

_pVS? 
M <*~ 4Iy 

_pSc 

Lq " Am Cl <* 

M - pVS ° 2 C 

r pSV r 

Ls ' 2 m Cl > 

„ pv 2 sa _ 

Ml 2Iy Cm ‘ 

C -*** 
a 2m a 

C -f> V2S ~ C c 

C "~ 2Iy q 

d -PSV n , 
"a - 2 m aa<>i 

R _ pV 2 S~c 

B< >- 2 Iy a " b ‘ 


Table II. - Characteristics of an advanced fighter aircraft 

and flight conditions. 

c = 3.51 m 

Cl 0 =2.7 

S = 37.16 m 2 

Cl* =2.5 

m - 15000 kg 

C Lq = 36. 

/y = 170000 kg-m 2 

C Lj =0.83 

p = 0.56 kg/m 3 

C ma = -0. 18 

V = 90 m / sec 

C m . =-2.5 


C m , = -10. 


C mj = -0.88 
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Table III. Estimated parameters from simulated data. 
Advanced fighter aircraft. 


Parameter 

True 

value 

Estimates ( ft ) 

No unsteady 
effects 

bj known 

bj estimated 
(c) 

c a 

2.75 


2.4 

2.75 




(0.23) 


a a 

-0.05 


0.4 





(0.25) 


bi 

1.0 


1.0 







Cl « 

36.0 

36. 

49. 




(3.7) 

(9.6) 

(1.2) 

C L S 

0.83 

0.7 

0.8 

0.829 



(0.10) 

(0.13) 

(0.0091) 

C L a 

2.70 

2.71 

2.8 0>) 

2.71 0 >) 



(0.043) 

(0.25) 

(0.062) 

MM 

84.35 

0.0828 

83.57 

0.0829 

83.61 

0.0828 

99.95 

0.0000172 


a Numbers in parenthesis are standard errors, 
b Computed from Cr =a a +c a . 

c Measurement noise reduced to 1/20 of its previous value. 
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Table IV. Estimated parameters from simulated data. 
Advanced fighter aircraft. 


Parameter 

True 

value 

Estimates (*) 

No unsteady 
effects 

bj known 

bj estimated 
(c) 

C <7 

-0.23 


-0.17 

-0.23 ! 





(0.064) 

a Q 

0.05 



0.04 




(0.038) 

(0.060) 

bi 

1.0 


1.0 

1 . 





(1.3) 

C m g 

-10.0 

-11.5 

-12. 

-10. 



(0.56) 

(1.4) 

(1.2) 

C m s 

-0.88 

-0.88 

k mm 

-0.879 



(0.016) 

mgsm 

(0.0091) 

Cm a 

-0.18 

-0.186 

-0.19 (M 

-0.2 (b) 



(0.0064) 


(0.12) 

R 2 (%) 

85.45 

83.00 

83.00 

97.82 

°{C m ) 

0.00124 

0.00124 

0.0125 

0.0000172 


a Numbers in parenthesis are standard errors, 
b Computed from C ma =a q +c q . 

c Measurement noise reduced to 1/3 of its previous value. 
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Table V. Correlation matrix of parameters (bj = 1.0) 


c a 


C L S 

a„ 

1.000 

-0.929 

-0.586 

-0.982 


1.000 

0.689 

0.922 



1.000 

0.565 




1.000 


Table VI. Correlation matrix of parameters 


c a 

a a 

bl 

Cl * 

C b s 

1.000 

-0.999 

-0.961 

-0.970 

- 0.778 


1.000 

0.952 

0.975 

0.783 



1.000 

0.879 

0.671 




1.000 

0.849 





1.000 


Table VII. - Characteristics of an experimental aircraft 
and flight conditions. 


c = 3.76 m 
S = 21.02 m 2 
m = 6700 kg 

1 Y =46792 kg-m 2 
q = 2211 Pa 
V = 101 m / sec 


C* La =1.12 
C L . =66.8 

< - io 

C*l 6 = 122 
C^ o =-0.142 
C TO . =-9.08 
C m<r =-1.29 

C* ms =-0.71 
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Table VIII. Estimated parameters from simulated data. 
Experimental aircraft. 


Parameter 

True 

value 

Estimates (®) 

No unsteady 
effects 

bj known 

bj estimated 

* 

2.55 


2.4 

2.8 




(0.11) 

(0.54) 

a a 

-1.43 


-1.2 

-1.5 




(0.14) 

(0.53) 

bj 

1.15 


1.15 

1.4 





(0.34) 

% 

1.0 

44. 

8. 

1 . 



(2.0) 

(4.7) 

(11.) 


1.22 

1.72 

1.25 

1.2 



(0.085) 

(0.099) 

(0.14) 

C l 

1.12 

1.50 

1.2(b) 

1 . (b) 



(0.025) 

(0.21) 

(1.1) 

R 2 {%) 

84.88 

82.03 

84.88 

83.25 

s { c l) 

0.0632 

0.0654 

0.0632 

0.0632 


a Numbers in parenthesis are standard errors, 
b Computed from Cjr, a = a a +c„. 
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Table IX. Estimated parameters from simulated data. 
Experimental aircraft. 


Parameter 

True 

value 

Estimates (®) 

No unsteady 
effects 

bj known 

6j estimated 

* 

c q 

-0.49 


-0.42 

-0.9 





(0.67) 

a q 

0.35 



0.8 





(0.66) 

bi 

2.07 


2.07 

3. 





(1.3) 

Cm q 

-1.29 

-9.4 

-3. 

2. 


! 

(0.34) 

(1.3) 

(6.7) 

Cm 6 

-0.71 

-0.79 

-0.72 

-0.69 


(0.014) 

(0.020) 

(0.038) 

c* 

^m a 

-0.14 

-0.166 

-0.1 (b) 

-0.1 <b) 



(0.0042) 

(0.11) 

(1.1) 

R 2 {%) 

s{C m ) 

85.31 

0.0109 

82.59 

0.0110 

83.00 

0.0109 

83.03 

0.0109 


a Numbers in parenthesis are standard errors, 
b Computed from C m ^ = a q +c q . 
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Figure 3. Computed time histories of aerodynamic coefficients with and 
without unsteady aerodynamic terms. Advanced fighter aircraft. 
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t, sec 


and predicted aerodynamic 
X-31A aircraft. 







t, sec 

Figure 11. Time histories of measured and predicted aerodynamic 
coefficients and residuals. X-31 drop model. 
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